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SKYLAB DATA AS AN AID TO RESOURCE MANAGEMENT 
IN NORTHERN CALIFORNIA (EPN NO. hSk, TASK ^J.3) 


.0 INTRODUCTION 

In order to determine the information content of Skylab S190 infrared 
ektachrome (CIR) imagery for wildland resource inventories, a timber inventory 
based entirely on manual photo Interpretation techniques was performed. In 
this study the SI 90A Imagery was utilized as the first stage in a multistage 
sampling design. The objectives included: 1) determining the abilities of 

human photo interpreters to distinguish timber volume classes on SI 90 data: 

2 ) determining the efficiency of simple random sampling applied to a multi- 
stage sampling design and 3) demonstrating a method of obtaining timber 
volume information of a reasonably large area (192 square miles located within 
the Plumas National Forest) that is timely, cost-effective and does not require 
the use of a computer. 

The need for timber volume surveys is based on the premise that sound forest 

management practices are dependent upon the availability of accurate, timely, 
and economical forest inventory data, which includes timber volume information. 
Forest inventory requirements in the Intensively managed temperate regions of 
the world include data not only on timber volume but also on timber stand 
condition, growth rate, ownership, soils, bedrock geology, surface geology, mineral 
extraction, subsurface water, surface water, vegetation, wildlife, land use, 
land productivity, climate, historical and cul tural - patterns , population market 
values, and transportation. However, 28% of the world's land area Is covered 
by forest, containing some 12^ trillion cubic feet of timber, and only a small 
portion of this vast area is intensively managed. Within the non-intensively 
managed or unmanaged forests, low-cost "volume only" surveys are needed, as 
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man for the first time finds it possible to place these forests under management. 
This project will demonstrate a procedure for providing low-cost timely "volume 
only" information without requiring any- sophi sticated computer terminals. The 
procedures for selecting sampling sites and estimating timber volCimes from the 
Skylab S190/\imagery and aircraft imagery at each stage of the sampling scheme 

will be discussed in following sections and will indicate the roles played by the 

S19QA imagery , supporting high-flight imagery, and large-scale aerial photography 
in performing this "volume only" survey. 

Gross timber volume was the variable estimated in the sampling design based 
on random sampling at each of three stages. Sampl ing. uni ts were selected at 
each stage at random because it was not known how well human interpreters could 
differentiate between timber and non-timber classes on S190Aimagery in the first 
stage, and by utilizing random sampling techniques, the ability to separate 
these classes could be established quantitatively. In the first stage timber 
'volume was estimated from the percentage of area occupied by each of 4 timber 

volume classes within 2 mi . x 2 mi . squares delineated on the S190A imagery . For 

the second stage, randomly selected squares from the first stage were transferred 
by visual inspection on to 1:60,000 and 1:120,000 infrared ektachrome photography 
and then divided into five rectangular flight lines. Photo interpreter delineations c 

the timber volume classes were measured within each square and areal percentages 

of the volume classes were calculated. A number of these rectangular blocks 

were selected at random for the acquisition of large-scale high-resolution 

aerial photography coverage. Measurements of percentage crown closure and 

average stand heights were made on .k acre plots located on the large-scale 

photography to arrive at the third stage photo volume estimate using Chapman's 

photo stand volume equations. Individual photo plot volumes were calculated 

-'Chapman, Roger C. 1965. Preliminary Aerial Photo Stand Volume Tables for 
Some California Timber Types. USDA-FS, Pacific Southwest Forest and Range 
Experiment Station, Research Note 93* 
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from regression estimators derived from precise measurements of volume made on 
randomly selected' trees with a Barr & Stroud dendrometer during the ground 
phase of the study. The ground volume measurements were then expanded through 
the various stages of the sample design to estimate the total timber volume 
located on the 192 square mile study area (see Fig. 1 for a diagram of the 
sampl ing design) . 

2.0 PROCEDURE 

2. 1 Selection of a Study Area 

An area, closely approximating that found on the USGS Bucks Lake 15 minute 
quadrangle, was located on a Skylab 3 SigOA infrared ektachrome transparency and 
enlarged to a scale of 1:170,000 on an 8 x 10 inch print (Fig. 2). The area, 12 
miles by 16 miles on the ground, was chosen because of historical and current 
information pertaining to it that the Remote Sensing Research Program currently has 
on file. Also, complete coverage of the area with recent 1:60,000 and 1:120,000 
resource photography was available. The SISOA infrared ektachrome imagery was 
chosen for interpretation because the contrast betv/een timber and non-timber areas 
is more readily apparent especially hardwoods and brush versus conifers. 

2 . 2 Volume Estimates by Human P.I.'s on Primary Sample Units (PSU's) 

The 12 mile x 16 mile study area was divided into ^8 squares (primary 
sample units: PSU's), each 2 miles on a side. The size of the PSU's was based 
on (l) a practical area that can be interpreted on S190A imagery enlarged to a 
scale of 1:170,000, (2) an area easily divided into flight lines which can be 
photographed by a. light aircraft using a 35mm camera system and (3) the ability 
of the ground crew to complete the ground work for a flight line in one day. 

Within each square, interpreters estimated the percentages of the area occupied 
by, 1) a timber density of less than 10,000 bd.ft./ac = T1 , 2) a timber density 
of 10,000 - 20,000 bd.ft./ac = T2, 3) a timber density greater than 20,000 bd.ft./ 



_/ 



FIG. 1 .DIAGRAM 'OF THE . 


A. A portion of the enlarged Skylab 3 
SI 90A I nf rared Ektachrome image 
showing the PSU grid for volume 
estimates for Stage I. 


B, A selected PSU transferred to a 
1:60,000 resource photograph and 
divided into 5 rectangular flight 
lines for the volume estimates for 
Stage II to be made on each flight 
1 ine. 


C. Ten .4 acre circular photo plots 
located on a selected flight line 
for the Stage III volume estimate. 


D. Selected photo plot to be visited 
on the ground for precise tree 
volume measurements. Encircled 
tree is one of four selected to be 
measured for volume data input to 
the regression estimator formula 
that will- estimate individual tree 
volumes on all photo plots. 


ING DESIGN 
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Fig, 2. The above image shows the 12 miles x 16 miles study area enlarged 
to a scale of 1:170,000 from a Skylab S190A Infrared ektachrome photograph. 
The grid represents the 2 miles x 2 miles Primary Sampling Units. 
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ac, = T3, and h) non-timber = NT (brush, hardwoods, bare soil, rock, grass, 
etc.). The Interpreters based their estimates on comparisons with training 
areas chosen from 1:120,000 scale color infrared photography acquired by the 
RB 57 aircraft. 


each volume class was multiplied by a specified weighting factor. For a 
volume of less than 10,000 bd.ft. per acre, the weight was 1. The factor for 
10,000-20,000 bd. ft. per acre was 2 and for greater than 20,000 bd.ft. per acre 
the weight was 3. Non-timber classes were multiplied by 0. The timber-class 
weighting factors were chosen as representative of the relative magnitude of 
volume present i.e. multiplying the area of a PSU in TT by 1 gives a relative 
volume figure for the area in class T1 . 

2 . 3 Selection of Primary Sample Units (PSU's) for Secondary Sample Unit (SSU) 

Volume Estimatio n 

The next step was the selection of a number of the PSU's for partial 
coverage with large-scale high resolution aerial photography and the delineation 
of second stage volume estimates. Six squares were chosen at random from the 
popuiatio:'. of A8 PSU's and delineated on high altitude 1:60,000 photography 
(Fig. 3a). The sample size was based on the variability of the area and upon the 
uncertainty as to how well volume classes can be interpreted on Skylab SI 9 OA 
infrared ektachrome imagery. 

Each of the six squares was divided into five flight lines. Each 
flight line was a rectanguiar area 2 miies long and approximately 2,100 feet 
wide, or iented in an east -west direction (Fig. 3b). One flight line was 
selected at random from the five possibie in each square. The volume 
estimate for this flight line was made in a fashion similar to the volume 


To arrive at a volume estimate within each PSU, the percentage area of 
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Fig. 3a. The image above illustrates one of the PSU's (PSU No. 3) delineated 
on a portion of a 1:60,000 high altitude color infrared photography. The 
boundaries of the PSU have been carefully transferred from the S190A enlarge- 
ment and, due to relief displacement, no longer appear square. 
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Fig. 3b. This image shows the same PSD with the four timber volume class 
delineations present as well as the divisions for the secondary sampling 
units. The dotted line in the center of one of the SSU's indicates a 
randomly selected flight line for acquisition of low-altitude high resolu- 
tion photography for tertiary sample unit volume estimation. 
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estimate for the PSU's on the S190 imagery. The only change, however, 
that the volume classes were delineated within each block and a coordinate 
digitizing device connected to a mini-computer was used to calculate the 
area of each class present. From these area calculations, percentages 
of total area occupied by each volume class within each chosen flight 
line could be determined and the appropriate weights applied for calcu- 
lating a total volume estimate for each selected flight line (SSU) . 

Each selected flight lliie was also located on a topographic map sheet 
to determine its average elevation above sea level. This information 
along with the flight line del i neatlon on a high altitude 1:60,000 photograph 
was used as navigational aids In acquiring the large-scale high- 
resolution aerial photography used to make the third. stage volume estimations. 
2.k Lar ge-Scale High Resolution Aerial Photography Acquisition and 
Interpretation 

Two 35 mm cameras were used to obtain low altitude photography of 
the selected secondary sampling units simultaneously at two different 
scales (Fig. ^) . A 24 mm focal length, wide-angle lens was used to ac- 
'•omplete coverage of each sampling unit at an approximate scale of 
1:7,500, and a 200 mm focal length was used to obtain large scale stereo 
triplets, scale approximately 1:1,000, from which to make precise photo 
esti.mates of timber volume. The camera with the telephoto lens vias equipped 
with a motorized film drive which enabled each stero triplet to be taken 
within one second at five second intervals while the camera with the wide- 
lens was operated manually to obtain single frames at five second inter- 
vals. The photo coverage for each SSU consisted of ten stereo triplets 
and ten wide angle photographs. 
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Fig. h. The top photo was taken with a 2k mm lens and the bottom photos 
illustrate the stereo triplet obtained by the telephoto lens at the same 
time. The inked circle on 'the bottom center photo delineates the scale- 
adjusted 0.4 acre ground plot. 
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The wide-angle photos of each secondary sampling unit were mosaicked 
together to show its full area. The center of the middle photo for each 
stereo triplet was used as the plot center, and these were located and 
marked on the mosaic. The plot centers were also located on a topographic 
map and the elevation of each was determined. 

The scale of each photo plot was determined, and a O.^-acre circular 
plot was drawn about the photo plot center. For all photo plots, average 
stand height and percentage crown closure were estimated in order to make 
individual photo plot volume estimates (calculated by Chapman's stand photo 
volume equations) that serve as the tertiary sample unit (TSU) estimate. 

Two plots were drawn at random from each SSU to be visited on the ground for 
precise plot volume measurements, i ncl udi ng four trees selected at random 
from each plot for measurement with an optical dendrometer. 

2 .5 Ground Work 

Each tree selected above was measured precisely on the ground with 

a Barr & Stroud optical dendrometer to calculate the individual volumes 

(Fig. 5).* The large scale photographs were used to locate the photo plot 

centers as well as the trees within the plots to be measured. All trees 

within the plot were identified as to species and measured for DBH. The 

plot radius as located on the ground was measured for its true distance 

and this was used to more accurately estimate both the photo scale and 

the plot area. The dendrometer measurements were brought back from the 

field and entered into a computer program that calculated merchantable 

stem volumes by species for the individual trees. These volumes were 

*The measurements taken with this device consisted of diameters outside bank, 
at specified points along the stem of the tree, in order to determine tree 
vol ume. 
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Fig. 5. The photo on the left shows the Barr and Stroud optical dendro- 
meter in use. The photo on the right shews the tree diameter being 
measured with a diameter-tape. 
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then regressed with photo volume estimates to derive the prediction equations 
for stand volume given a photo volume estimate. The volume estimates derived 
at each stage were expanded through the sample design (multirstage random 
sampling) to estimate the total volume for the study area. 


2 


I 


3.0 RESULTS 

Since selection of sample units was With equal probabilities, 
nn.;qible t^ n<;inn two different estimators and also using__t_wo 

d ifferent auxiliary var_Lab,le.S-.. Using a ratio estimator allowed combining 
the information available from the SlSOAimagery with direct measurements, and 
for comparative purposes a simple random sampling estimator was also used. 

The auxiliary data consisted of a'weighted volume estimate, described 
earlier and also the number of acres In timber, obtained from the basic data 
described earlier. Results of these estimates are portrayed in Table 1. 

While the estimates obtained by each of these methods are roughly comparable 
in magnitude, the most important observation is that simple random sampling 
was actually more precise,with an .RSE of 0.1315'. - This result may at f.irst 
seem strange since this technique does not take advantage of the auxiliary 
SlSOAdata; -however as is pointed out by Cochran* auxiliary variables do not 
necessarily lead to increased precision. Gains in precision depend on the 
correlation between the auxiliary variable, x, and the variable of interest, 
y. Cochran gives a formula for a minimum correlation, Cq, which must be 
met in order for gains in precision to occur. By comparing the correlations 
obtained' in this application with Oq in Table 1 it is easily confirmed that 
no gain in precision is likely to occur with this particular auxiliary data. 
Therefore, for this particular population and S190A "Imagery, i t- can be 
concluded that simple random is the most precise among these three estimation 

/ 

alternatives. Details of the estimation procedure are given in Appendix 1 . . 

' •' 

f 


*Cochran, W.O. 1963, Samplmg Techniques, pg.;l65. 
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Table 1 Summary of Estimates by THREE Estimators 


^\Estimatot 

Estimates^^^ 

Ratio 

X = Wtd . Vol • Est . 

Ratio 

' X = Area of Timber 

Simple 

Random 

Y 

528,663,956 

517,386,458 

547,838,296 

se(Y) 

83,108,888 

.81,256,237 

72,054,597 " 

RSE 

0.1572 

0.1571 . 

0.1315 

✓V 

R 

2891 .9 

5802.6 


Pxy 
R S 

-0.1766 

0,0395 

- 

p — ^ ' 

2 S 

y 

1 

0.2124 

0.2811 



DEFINITIONS: 


Y 


= Estimated total volume of timber In cubic feet trees 5" DBH 


SE(Y) 

RSE 


= Standard error of Y 


= RELATIVE STANDARD ERROR, SE(Y)/Y 


^ A. 



R 


= 2y/2x 


p = Correlation between x and y 

xy 




and S = Standard deviation of x and y, respectively. 




Skylab S130B Imagery was not available for this study area so It was not 
possible to make any direct comparisons between the S190A and SI 9 OB products. 
Inspection of the higher resolution SI90B imagery covering an analogous 
area, however, indicates potential gains in inventory precision if it is used 
In place of the S190A imagery in the weighted estimate multistage design. 
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Estlmation Procedure 

As described earlier, the sampling system Is comprised of three stages with 
random sampling, without replacement, at each stage. Since the population is 
finite and relatively small, the finite population correction is appropriate. 
Following the notation of Cochran 1363, Sampling Techniques (pg. 285 ) ;the values 
for y, the estimated population mean volume per acre, and its estimated variance v(y) 

are obta i ned 'aS‘ fol lows . 


Let 

t h 

y = volume/acre as "measured" at the ijk .4 acre circular plot 
^ i jk 

.th ,th . th 


v~7V f CChapman volume estimate) = Vol . for the I PSU, j SSU, k 

TSU obtained by regressing sample data from both photo and ground 


- . . data- 
Then, 

= n m k 

y = 222 y.../ nmk , 


i j k 


ijk 


And 


y 


v(y) = 


ilA 2 

n ^ 1 


6 = n first, stage samples 
1 = m second stage samples 
10 = k third stage samples 

y 


- V 2 * , - ^) 2 

nra ^ 2 ■ nmk , ^3 


^ - £L-2_ - IOC f _«l-l- o 

^1 - N - 48 > ^2 M “ 5 “ 


' h 


~ = .0015625 

(Sampling fractions) 


2i.v,-y 2 


n - 1 


n - 1 







I- I, 




V \ 

) 




^ ,u, K-r 'c-'C.-^^y' . 
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n m 


I j IJ I I J 'J ' 


n m 




-2 = 2 ^ 


n(m - i) 


n(m - 1) 


Y; = 


m k y/ 

^ 2 y. . , 

•, "ijl 


mk 


n m k 

2 2 2 (y.j^ - y.j) 


n m 


i J 1 


nm(k - 1) 


2 2 ( 2 Y- . 1 - ky . . ) 

• i J ^ •) U ^ 

nm ( k - 1 ) 


"^ij 




* Not available since m - 1 . Hence estimated variance is lower than 
true variance. Justification can be given for using the approximation with- 

i 

out second term since f^ = .125 and in other applications f^ might be even 
smaller. Thus the first term asymptotically represents the primary con- 
tribution to the sampling variance, and the effects of later stages of 
sampling are reflected in the value of the estimate rather than its form 
(Kendall & Stuart, The Advanced theory of Statistics , Vol . 3, P- 200). 


Since y estimates the average volume per acre, the total volume 


on the area (Ty) may be estimated as 

Ty = Xy , where X = acreag e of populatio n— 

A 2 - 

v(7y) = X v{y) , assuming that the acreage is known wi thout- error. 


Confidence intervals are calculated in the usual way, 

\/ A ■ A \f~A 

Vuifil < ft < fl +.1-.-^ . .. Vv(0 


- ^^-1’ a/2 


ct/2 


(0) ) - 1 - a 


or 


A V / 

6 + t Vv(0) with probability level 1 - a 

— m-K a/2 


where 


0 = the population parameter of interest 

A - A ' 

0 “ the estimate of Interest (either y or Ty) 

A. ^ A 

v(0) “ estimated variance of 0 
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APPENDtX I! ANALYSIS OF THE INFORMATION CONTENT OF SKYLAB S190A DATA 
Correlation Analysis as Applied to Photo-Interpretation Results 

' The first stage in this multi-stage inventory design v;as for a human 
interpreter to separate three timber volume classes from non-timber on en- 
larged SI 90 A imagery. In order to determine how accurately this could be 
done, 18 of the 48 total PSU's were delineated and interpreted for volume 
classes on Zeiss scale 1:60,000 infrared Ektachrome high altitude imagery. 

The 18 selected were the only ones for which stereo coverage was already 
on hand. A skilled interpreter delineated the 3 timber volume classes and 
the non-timber class on the high altitude imagery and precisely determined 
the percentage of the total area covered by each class with an ..aj:ea--caI.CH- 
laf inn device. These results, assumed to be an accurate account of ground 
conditions, were correlated with the interpretation results of each of three 
interpreters on the enlarged portion of the SI 90A' imagery. “The initial 
analysis (Tables 1, 2, S3) showed poor correlation (generally less than .50) 
between interpreters and high flight imagery over all volume classes,* 
including the non-timber class, s.uacLe^tJLnq Lt_Ls_Jlo.t -P-as.sJLbJLe to differen- 
tja_|e,,volu^,_,cl^^_ses_on__the^J^ volume classes of 10,000- 

20,000 board feet/acre and greater than 20,000 board feet/acre were combined in 
an attempt to improve the correlation but the results showed improvement only 
In the non-timber class for all interpreters (.62-. 67 correlation) while the 
timber volume classes ranged from .44 to -.07. This only strengthened the 
conclusion of the Inability to differentiate timber volume classes. However, 
when all volume classes were combined to see the correlation between timber 
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and non-timber as Interpreted onSlSOA and as delineated on 1:60,000 photo-, 
graphy (Table A) the coefficients ranged from .62 to .67 for non-timber and 
from .61 to .66 for timber. The correlations between interpreters on SI90 
imagery alone (Table 5) were .81 to .91 which leads to the conclusion that it is 
possible to consistently differentiate between timber and non-timber on 
enlarged SI90A imagery , but it is not possible to consistently differentiate 
any volume classes within the timber -areas . 

Analysis of Variance and Multiple Comparison Tests 
To Judge Significance of Photo- Interpretation Results 

In order to quantify and determine the significance of error for indivi- 
dual photo- interpreter (P.l.) results among timber strata, a-_ two-way analysis 
of variance (Scheffe, 1359) was performed. The three Skylab photo- interpreters 
(P.I.'s) were considered levels of factor two. Data for each P.l. -timber 
stratum combination consisted of the difference between the appropriate Skylab 
and high-flight percent timber area estimates on .each of the aforementioned 
18 PSU‘s. wag agc;tifnp.rl tha t_tba.s e I fff inHiviHual PSU ' S Were 

iJLde,p_ejidej3JiJ^LtJ3lD_aj3^between_p.LcLtPj^^ overall distri- 

bution Q.^the_dl££er.en.ce3-.aP - PX.oaclie.d_a . n.Qm distribution. 

The analysis of variance then gave least-square .e.5-t,i mates for the 

^CCOJI. (relative to the high-flight data) due-JLo using ( 1 )_ a given-photo- 
interpreter over all timber strata, (2)_ a given timber stratum over all 
P.I.'s, and a given P.i. -timber stratum combination. These error esti- 
mates are given in Table 6. F-statistics were then formulated (see Table 7) 
from the variance of these expected errors to test the null hy potheses that 
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TABLE 1. CORRELATIONS^ BETWEEN PHOTO- INTERPRETER A ON SKYLAB SI 9 OA IMAGERY 

AND PHOTO- INTERPRETER D ON ZEISS 1:60,000 HIGH-FLIGHT PHOTOGRAPHY 
BY TIMBER STRATUM 


Photo- I nterpreter A 



- 

T1 

T2 

1 

T3 

T4 




^T1 \ 


/. 69/'' 

-.30 

-.il6 

T1 

= Timber Stratum 1 

Photo- 

■ T2 

.25 


-.05 

-.23 

T2 

= Timber Stratum 2 

1 nterpreter 








D _ 

r 

T3 

.01 

-.26 

<<3 

. .07 

T3 

= Timber Stratum 3 


NT 

-.21 

-.67 

.27 

(3^ 

NT 

= Non-Timber Strati 

TABLE 2. CORRELATIONS 

BETWEEN 

PHOTO- II 

UTERPRETER B ON 

SKYLAB : 

SI 90 

IMAGERY AND 


PHOTO- INTERPRETER D ON ZEISS 1:60,000 HIGH^FLIGHT PHOTOGRAPHY, BY 
TIMBER STRATUM 


Photo- I nterpreter B 




T1 

T2- 

T3 

't4 




•T1 

.33 

.39 

-.14 

-.64 

T1 

= Timber Stratum 1 

Photo- 

T2 

.16 

-.12 

-.24 . • 

.06 

T2 

= Timber Stratum 2 

I nterpreter 
D 

■*■3 

-.35 

.09 

:i2 

.15 

T3 

= Timber Stratum 3 


NT 

-.26 

-.35 

.35 

.49 

NT 

= Non-Timber Strati 

TABLE 3 . CORRELATIONS 

BETWEEN 

PHOTO- j 

INTERPRETER C ON 

SKYLAB 

SI 90 

IMAGERY AND 


PHOTO- I NTERPRETER D ON ZEISS 1:60,000 HIGH-FLIGHT PHOTOGRAPHY, BY 
TIMBER STRATUM 


Photo- Interpreter C 




T1 

T2 

T3 

t4 



T1 

.49 

.54 

-.20 

.66 

T1 

Photo- 

T2 

.03 

-.05 

-.10 

,.19 

‘'12 

Interpreter 




- 

» * ' 


D 

T3 

-.22 

-31 

-.28 

.02 

T3 


NT 

-.41 

-.75 

.54 

.47 

NT 




Timber Stratum 1 
Timber Stratum 2 
Timber Stratum 3 
Non-Timber Stratil 
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TABLE h. 


CORRELATIONS BETWEEN PHOTO- INTERPRETERS ON SKYLAB S190A IMAGERY 
AND PHOTO- I NTERPRETER D ON ZEISS 1:60,000 HIGH-FLIGHT PHOTOGRAPHY 
CONSIDERING ONLY TIMBER VERSUS NON-TIMBER INTERPRETATION 


Photo- 
I nterpreter 
D 


Photo- 1 nterpreter A 

Timber Non-Timber 

Timber ■ .66 -.67 

Non-Timber -.66 *67 



Photo- 1 nterpreter B 



Timber 

Non-Timber 

Photo- 

Timber -61 

- .62 

nterpreter 

D 

Non-Timber *’•61 

.62 


Photo- 
i nterpreter 
0 


Photo- 1 nterpreter C 
Timber 

Timber . 63 ' 


Non-Timber 

.64 


Non-Timber 


-.63 


.64 



TABLE 5. CORRELATIONS 
CONSIDERING 

BETWEEN PHOTO- INTERPRETERS ON 
ONLY TIMBER VERSUS NON-TIMBER 1 

SKYLAB SI 30A IMAGERY 
INTERPRETATION 


Photo- 1 nterpreter B 



Timber 

Non-Timber 

Photo- 

Timber *81 

-.81 

Interpreter 

A 

Non-Timber -.81 

.81 


Photo- I nterpreter C 



Timber 

Non-Timber 

Photo- 

Timber .84 

1 

oo 

I nterpreter 
A 

Non-TImfaer -.84 

oo 


Photo-Interpreter C 



Timber 

Non-Timber 

Photo- 

Timber -91 

-.91 

interpreter 

B 

Non-Timber “-91 

.91 
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TABLE 6. ESTIMATED EFFECTS (ERRORS) FOR FACTORS AND INTERACTIONS IN THE 

TWO-WAY ANAirY'SIS OF VARIANCE FOR SKYLAB PHOTO- INTERPRETATION RESULTS. 


FACTOR OR level LEAST-SQUARE ESTIMATE OF 

INTERACTION FACTOR OR INTERACTION EFFECT 


Photo- 1 nterpreters 

P.I. A 

a 

1 

— 

-12.95 

- 

P . I . B 

a 

2 

= 

15.88 


p.r. c 

Ay 

a 

3 

z= 

- 2.93 

Timber Stratum 

1 

e 

1 

=r 

- 1.34 


2 

3 

2 

“ 

.012 


3 

/V 

3 

3 


1.328 

Photo-Interpreter x 
Timber Stratum 

-A-x 1 

A V ? 

/V 

V 

= 

- 3.07 
10.86 

- 7.79 


/% 

A X 3 

^2 

y\ 

y 

3 



B X i 

Y 

1 

= 

- 3.51 


B X 2 

^2 


- 9.48 


B X 3 

/% 

^3 


12.99 


Cxi 

> 

1— ■ 


6.58 


C X 2 

yv 

Y 

2 

= 

1.38 


C X 3 



- 5.20 


>-i “ « "• Pi . , 



- 2 ^»- 


TABLE 7 . F-TESTS BASED ON TWO-WAY ANALYSIS OF VARIANCE 


For P,l. effect (error) over all timber strata 


sum of squares - sum of squares ^ 1 2 

, 7 1 V = l '5 

83893 - 60750 . Jii 
60750 2 


F = 29.14 p <.01 


p in?. 


2. For timber stratum effect (error) over all P.l.‘s 

^ _ sum of squares ,g - sum of squares ^ ^ , 2 

Sum of squares^^^^^ V = 15 


60942 - 60750 . 153- 

60750 2 


F = .24 p >.10 

For Interactions effects (errors) over all P.l.'s and all timber strata 
sum of squares - sum of squares ^ ^ ^ ^ 

p = - - q V = 15: 

Sum of squares 


70300 - 60750 . iii 

60750 4 


6.02 


p<. 01 


I “i. 45 



^ P i-r» r g The P. 1 . and. interaction hypo 


theses v.»ere rejected at the_0CLeji[: 


Indicating that^sona^ 


timber strata hypothesis was accepted a, 


and i t was 


therefore concluded that Dhoto-interp_r.e-t-a±.i.o n. .ftrrnr_s_attr i buJbAble. .to differ - 
ftnr.fts in timber Strata v/ere not significant^ 

To determine which P.l. errors were significant t JitLe.e__ty.e.e.s_p f 
mi I If- i P lA-e:Qmpa,t;.i:^ n»fees=<?g (Scheffe, 1959) were applied to the results of 
the analysis of variance. These were (1) the Tukey test assuming inde- 
pendence of P.l. errors, (2) the Tukey test' not assuming independence of 
P.i. errors, and (3) the Scheffe multiple comparison method. Since the 
analysis of variance requires the sum of the level effects for a given 
factor (i.e., the sum of the P.l. expected errors)to add to zero In order 
to obtain unique least-square estimates for those errors, the s econd Tukey 
test was included to check the first for any changes resulting from the 
violation of error independence. 

Parameters to be tested consisted of the difference between each 
calculated P.l. expected error value and its ideal value of zero. In addi- 
tion, differences were formed pairwise between interpreter errors to deter- 
mine which P.i. results significantly differed. All three multiple compari- 
son tests allow a statement of statistical significance for a given error 
difference to hold true for all such differences of the same magnitude for 
any chosen level of significance.-: 

The results of these three tests are summarized In Tables 8 through 
10. Note that all three multiple comparison tests give the same results 


ORIGINAL PAGE-IS 
OP POOR QUALITY 


- 26 - 


TABLE 8. TUKEY MULTIPLE-COMPARISON TEST (ASSUMING EFFECT INDEPENDENCE) FOR 
SIGNIFICANCE OF. P.i. EFFECT (ERRORS) ESTIMATED FOR THE TWO-WAY 
ANALYSIS OF VARIANCE 


Contrast (i.e. 
Difference) 

Test 

Statistic^ 

Computed 
Level of 
Signif icance^ 
(Reject Null 
Hypothesis of 
No Significant 
Difference if 

p< . 1 o) 

/Si A yN. 

Y =a '*a 
112 

q - -10.37 

p<. 01 

AAA 

2 13 

q = - 3.59 

.01 p<.05 

AAA 

T =a “a 
3 2 3 

q = 

6.74 

p<. 01 

A A 

=a “0 
4 1 

q = - k.6k 

p<. Ol 

A A 

'f “a ”0 
5 2 

q = 

5.69 

p <. 01 

T =a "o 
6 3 

q = -1.05 

p>.10 


where: 




3,153 


. 05 


3,153 


. 01 


2.93 

3.36 

h^20 


(1) Tukey Test Statistic “ 9 . 
Assuming Independence ’ 


as 


where- T = estimated contrast 

Inverse of the no. of observations made by a given 
P.I. in the experimental layout for the given 
analysis of variance -procedure 


s = square root of the mean square for residual error 
in the analysis of variance 


*^p = significance level 

k = ho. of P. !.. error’ estima.tes 


V = degrees of freedom for the residual error in the 
analysis of variance 


(2) Level of significance = p= probability that the null hypothesis of no significant 

difference. Is true given the computed contrast. 
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TABLE 9. TUKEY MULTIPLE-COMPARISON TEST {NOT ASSUMING EFFECT INDEPENDENCE) FOR 
SIGNIFICANCE OF P.l. EFFECT (ERRORS) ESTIMATED FOR THE TWO-WAY 
ANALYSIS OF VARIANCE 


(1) T"'''"/ 



2.93 

3.36 

^.20 


where T, a, s, p, k, v are as defined in Table 8, and 

A. 

cov (a. ,a:-) 

b ! 

2 

for the i and j of the contrast' of Interest, ij^j 


(2) Defined in Table 8 



TABLE 10. SCHEFFE MULTIPLE-COMPARISON TEST FOR SIGNIFICANCE OF P.l. EFFECT 
(ERRORS) ESTIMATED FOR THE TWO-WAY ANALYSIS OF VARIANCE. 


Contrast (i.e. 
Difference) 

Test 

Statistic^ 

Computed 
Level of 
Significance^ 
(Reject Null 
Hypothesis of 
No Significant 
Difference if 

p<.10) 

A A 

^ =a -a 

F = 18.79 

j 

p<. 01 where: F =2.35 

112 


2,153 10 

A ^ A 

^ =a -a 

F = 3.^1 

.01 p<.05 F ‘ = 3.07 

2 13 



A A /S 

'f =a -a 

F = 12.03 

P<.01 F = 4.79 

3 2 3 



A A 

^ =a “0 

F = 11.41 

p<. 01 

4 1 

' 


A A 

¥ =a -o 

F = 17.15 

p<.01 

5 2 

* 


A A 

^ =a -0 

F = .58 

p> . 1 0 

6 3 




(l) S<~h“ff4 Test Statistic = F 


'V 


p;q,v 


q? 




where. ^ = estimated contrast 

= estimate of the variance for the es.tLmated contrast 
'T - 9 c, 


1 


'^-1 


where s^'= mean square of residual error in the 
analysis of variance 

1 = index for the estimated P.l. effects ■ 

= no. of observations made by a given P.'i. in 
the experimental layout for the given analyst 
of variance . ‘j 

= a contrast defined b'y- the contrast of interej 

^.=2c-a . ■ 

■ 111 


'1 


1 


9 = degrees of freedom for P.l. effects 
q = no. of P.lv effects minus one 

V = degrees' of freedom- for the'fesidual error in the 
analysis of variance • 
p =• level of significance'- 


C2) Defined In Table 8. 
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except for the difference between the P.l. A and C errors In the second 
Tukey method. In this case, the difference is just slightly more statis- 
tically different from zero than indicated by the other two methods. 

Examination of Tables 8 through 10 shows that ^ 

had errors sianif icaatl y, djjf_f_e^ -_ent f rp_m ,zfii;o_over st_ra.t.a. That 

is, their results significantly differed from those obtained by interpre- 
tation of high-flight imagery. Only P.l. C. had no s TaiiTf icant,djj££e.Cfiac.e. 
from zero (at the 10 percent leve l) Indicating that the- expected Skylab 
interpretation error, relative to the high-flight base, attributable to 
that P.l. was not statistically significant. Note also that Jces.u.lJ:5.^ - 
a ll nhntn- interpreters over aLL tJ^lxec...sJxaJLa_.siqn i dif£ecejj_Lr om 

eg£tuSi±ietG««- The strongest differeaces occurred between P.l.'s A and B, 
and B and C, 

Differences involving expected errors due to P, I, -timber stratum 
interactions were formulated in a way s imi lar to ;that for P.L^s alone- 
These differences were then tested for statistical significance by the Scheffe . 
method. Results are given in Table 11. Examination of the table indicates 
that no interaction errors b,v_ themselves, when compared to zero error — axe 
tci-ioai iy g-Lq.n-i-f i oar>< - . How ever, three differences between interactio n 
ftr pors are s i a n i flcaii t ♦ Two of these involve photo- interpreters A and B, 

.one. in the context of timber stratum 2 and the other in timber stratum 3- 
■ The.n^ysical significance to be attached to this result may be obtained by 
examining the magnitude of interactions. Then from Table 6 it is seen that 
P.i. A tended to call more of timber stratum 3 as timber stratum 2 while 
P.l. Z did the opposite. 
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TABLE 11. SCHEFFE MULTIPLE-COMPARISON TEST FOR SIGNIFICANCE OF INTERACTIONS 
EFFECT (ERRORS) ESTIMATED FOR THE TWO-WAY ANALYSIS OF VARIANCE. 


( 1 ) 


Contrast 

Test 

Statistic^ 

Computed 
Level of 
Significance^ 
(Reject Null 
Hypothesis of 
No Significant 
Difference if 

p<. 10) 

^ > 

II 

> 

1 

O 

F = .11 

p>.10 where; F =1.99! 

1 1,1 
A. 



=Y -0 

F = .14 

p>.10 F = 2.45 

2 1,2 



II 

o 

F = .49 

p>.10 F = 3-48 

3 1,3 



11 

1 

o 

F = 1.34 

p>. 1 0 

4 2,1 



o 

II 

F = 1.02 

p>. 1 0 

5 2,2 



11 

-<> 

1 

o 

F = .02 

p>* 1 0 

6 2,3 



0 

1 

<>- 

11 

F = .69 - 

p>. 1 0 

7 3,1 



o 

J 

<>■ 

II 

F = 1-91 

\ / p>* 1 0 

8 3,2 



$ =Y 

F = .31 

p>*10 

9 3,3 



y\ /s 

T =Y “Y 

F = 2.34 

. 05 10 

10 2,1 2,2' 



^ A A 

Y =y-Y 

F = 2.45 

.05<P <• 10 

11 



A A A. 

Y =Y-Y 

F = 2.86 

.05"^5<.10 

12 







ScheffS Test Statistic = ^ 

where = are as defined as in Table 10, and 

q = degrees of freedom for interactions 
= no. of interactions minus one 
1 = index for the estimated interaction effects 
J =■ no. of observations in. the experimental layout for 

the given analysis of variance relating to the calcu- 
lation of the estimated interaction of interest. 

= a constant defined by the contrast of interest, 

T =2c Y ■ 

111 . 


(2) Defined on Table 8 
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In order to gain a better perspective on the effect o^ error associate d 
with each Dhoto-Interpreter-jaLLthl D-each. Xlmb_e.tij/.oJLum.eL-5Jtl^^^^ two types of 
on ft~wav analysis of variance first consisted of 

determining the size of expected errors for P.l.’s within each timber stra- 
tum and also the non-timber stratum (see Table 12). Data for this case 
consisted of the difference between Skylab and high-fllght timber stratum 
area estimates for each P.l. in each of the 18 PSU's which had high-fllght 
timber estimates. F-tests formulated from sums-of-squa res associated with 
P.I., effects, the remaining sums-of-squa res not accounted for by P.l. effects, 
and the associated degrees of freedom were used to test the null hypothesis 
that P.l. errors were not significant. 

Rejection of this hypothesis occurred only in timber volume strata 
1 and 2. A Tukey multiple comparison test assuming independence of P.l. 
errors was then utilized in these two strata. The assumption of P.l. effect 
(error) independence was utilized here since results from the Tukey multiple 
comparison tests based on the two-way analysis of variance results were very 
similar. Results are given in Table 13. In both timber strata 1 and 2, P.l. 

B and C effects were statistically significantly different, while A and C were 
as well, in timber volume stratum 1. Importantly, note that the F-tests for 
timber stratum 3 and non-timber and the Tukey test for 1 and 2 indicate that 
no P.l. errors differed significantly from zero. 

These last results differ with those of the second type of one-way 
analysis performed. Data for this second approach consisted of the actual 
photo- interpreter percent area estimates by PSU for. each timber stratum, 
in this case, estimates for the high-fllght were considered in the analysis 
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TABLE ]2. ESTIMATED P.I. EFFECTS (ERRORS) BY TIMBER STRATUM FOR THE ONE-WAY 
ANALYSIS OF VARIANCE FOR SKYLAB PHOTO- INTERPRETATION RESULTS 


LEAST SQUARE ESTIMATE 
LEVEL OF FACTOR EFFECT 


Timber 

Stratum i 

P.I. 

A 

Ak 

a = 
1 

k.k 



P. I . 

B 

a = 
2 

n.6 



P. 1 . 

C 

/V 

a = 
3 

7-2 

Timber 

Stratum 2 

P. 1 . 

A 

a = 
1 

6.09 



P. 1. 

B 

/V 

a “ 
2 

8.29 



P.I. 

C 

A. 

a = 
3 

IA.38 

Timber 

Stratum 3 

P.I. 

A • 

a = 
1 

7.33 



P.I. 

B 

a — ^ 
2 

4.60 



P.I. 

C 

a = 
3 

2.67 

Timber 

Stratum h 

P. 1 . 

A 

a = 
1 

- .60 



P.I. 

B 

A 

a = 
2 

1.73 



P. 1 . 

C 

As. 

a = 
3^ 

1.13 
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table 13. tukey multiple-comparison test (assuming effect independence) for 

SIGNIFICANCE OF P.l. EFFECT (ERRORS) BY TIMBER STRATUM ESTIMATED FOR 
THE ONE-WAY ANALYSIS OF VARIANCE. 


Stratum 

Contrast (i.e. 
Difference) 

Test 

Statistic^ 

Computed 
Level of 
Significance^ 
(Reject Null 
Hypothesis of 
No Significant 
Difference if 
P<.10) 

Timber 1 

^ A 

Y = a - a 
112 

q = -2.99 

,05<P<. 10 


\rr — a - a 
2 1.3 

q = .52 

p>.10 


ys yv 

Aj; = a - a 
3 2 3 

q = 3.51 

.oi<p<.05 


A, Ax 

w z= a - 0 
4 1 

-q = - .82 

p>. 1 0 


Ax Ax 

= Ot " 0 

5 2 

q = :2.17 

p>. 1 0 


w = a - 0 
6 3 

q - -1.35 

p>. 10 

Timber 2 

Ax Ax A 

^ = a - a 

112 

q = .36 • 

■. ■•-P>.10 


Ax Ax AX 

(j/ = a - a 

2 13 

q = - 3.49 

.01 <p<. 05 


Ax Ax Ax 

^ = a " a 

8 2 3 

q = - 3.89 

.oi<p<. 05 


Ax Ax 

f “ a - 0 
4 1 

q = .99 

p>. 1 0 


$ = a - 0 

5 2 

q = 1.35 

p>. 1 0 


Ax Ax 

f =s a - 0 
6 3 

q= 2.34 

p>.10 

Timber 3 

II 

q =■ 2.56 

p>.10 


AX Ax 

'j» = a - a 
2 13 

q = 2.14 

p>.10 


Ax /X Ax 

Y = a - a 

q = - .41 

p>. 10 


3 2 3 



Ax Ax 

Y = a - 0 
4 1 

q = 1.57 

p>.l0 


^ = a - 0 

5 2 

q = .99 

p>. 1 0 


Ax Ax 

^ - a - 0 

q = .57 

p>.10 


where: 


3,42 


.10 


3,42 


. 05 


3,42 


. 01 


2.987 

3.^36 

4.271 


-3k- 


Timber k 

4\ ^ 

T = a - a 
1 1 2 

q = “ .62 

p>. 1 0 

' 

A* /V A 

y = a - a 
2 13 

q = . 1 ^ 

p>. 1 0 


A A 

y ” a - a 
3 2 3 

q = 

.76 

p>. 1 0 


$ " a - 0 

4 1 

q = 

.16 

p>. 10 


y ^ a - 0 
5 2 

q = 

.kS 

p>. 10 


A A 

T = a - 0 
6 . 3 

q = 

• 30 • 

o 

A 

CL 


(1) Test statistic defined In Table 8 

(2) Defined In Table 8 
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as a fourth level of the photo- Interpreter factor. Effects were calculated 
(see Table 1^) and F-tests then rejected the null hypothesis of non-signlfl- 
cance of P. I . effects In each timber volume stratum. The Tukey test assum- 
ing Independence was then applied to the calculated P.l. effects. 

Examination of Table 15 Indicates that significant differences between 
the Skylab photo- interpreters generally follow the results of the first 
one-way analysis of variance. However, significant differences now exist 
between P.l.'s A, B, and C versus D (high-flight) P.l. effects. This last 
result is In contrast to the lack of significant differences relative to 
high-flight results found from the first type of analysis of variance dis- 
cussed above. From this last comparison it can be concluded that the 
F.l. effect based on PSD by PSU error (first one-way analysis of variance 
procedure) and the expected P.l. effect based on over-all error (second one- 
way procedure) tend to give different results. 

A comparison of photo- interpreter correlation results (Tables 1, 2, 
and 3) and photo- interpreter error (effect) sizes ' (especial ly from the one- 
wav analysis of variance) as depicted in Tables 6, 12, and lA allows the 
following observation: two interpreters' results can vary up or down together 
resulting in high positive correlations, but their respective average timber 
class estimates can be significantly different. This result suggests two 
d esixable objectives for photo- interpretation performance wPthin a given 
timbjer stratum from the standpoint of cost-effective sampling for timber . 
volume estimation. 

First there should be high .correlation between Skylab and high-flight 
estimates. This situation would a1 low meaningful stratification between 
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table ]k ESTIMATED. P. I . EFFECTS (ERRORS) BY TIMBER STRATUM FOR THE 
ONE WAY ANALYSIS 'OF VARIANCE FOR BOTH SKYLAB PHOTO 
INTERPRETATION AND HIGH-FLIGHT PHOTO . I NTERPRETAT I ON RESULTS 


LEAST SQUARE ESTIMATE 
LEVEL OF FACTOR EFFECT 


Timber Stratum 1 

P.l 

. A 

“i = 

6.08 


P.l 

. B 

/I 

= 

9.92 


P.l 

. C 

“3 

12.75 


P.l 

. D 

a/j = 

8.91 

Timber Stratum 2 

. P.l 

. A 

a, = 

1 .95 


P.l 

: ® 

“2 = 

^.95 


P.l 

1 . C 

“3 " 

18.52 


P.l 

1 . D 


1 1 .62 

Timber Stratum 3 

P.l 

1 . A 

/V 

8.82 


P. 1 

1 . B 

ys 

«2 = 

6.52 


P. 

I . C 

“3 

kA5 


P. 

1 . D 

“A " 

2.15 

Non-Timber Stratum 

P. 

! . A 

“l = 

-.75 


D 

i « 

1 . B 

a^ = 

1 .58 


P. 

1 . C 

A, 

“3 

1.35 


P. 

1 . D • 

“ 4 " 

.52 
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If' 

timber volume classes and would thus reduce the variance between PSU's and 
SSU's within a stratum. The result would be that fewe r PSU's and SSU's 
would need to be subsampled to obtain timber volume estimates within speci- 
fied allowable error and confidence level constraints. Therefore, costs 


of inventory would be minimized. ^ 

The second objective for photo- interpreter performance is Vow average 
expected errors for given timber strata. In this case, the probability 
that sample allocation will be to PSU's and SSU's that actually fal'1 in 
their appropriate strata will be maximized. The efficiency of the sample 
design in producing sample means for timber strata as close as possible to 
their true population means will then be maximized. 



